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Международный научно-исследовательский проект МиБоду Адуапсеа Айз$р Гог Ткапзрог (МААТ), финанси- 
руемый Евросоюзом по Седьмой рамочной программе, решает проблему планетарного масштаба — экологи- 
чески безопасной транспортировки. Один из разделов проекта предусматривает исследования аэродинамики 
специализированных высотных аэростатов — так называемых челноков. Они предназначены для обеспече- 
ния коммуникаций между наземными аэропортами и аэростатами, постоянно курсирующими в стратосфе- 
ре, — крейсерами. Решается одна из задач аэродинамических исследований, необходимых для получения 
соответствующей математической модели. Поставлена и решена задача математического описания закона 
изменения силы всплывания при вертикальном движении челнока в невозмущённом и возмущённом его по- 
ложении. Для этого проведён имитационный эксперимент, основывающийся на так называемой «электрон- 
ной обдувке» путём компьютерного моделирования. Методами регрессионного анализа обработаны резуль- 
таты компьютерных измерений силы сопротивления всплыванию в разных условиях. Полученное выражение 
позволяет идентифицировать математическую модель динамики движения челнока, представленную урав- 
нениями Эйлера. Результаты проиллюстрированы табличными и графическими данными. 

Ключевые слова: стратосфера, летательный аппарат, аэростатный принцип всплывания, силы и моменты 
сопротивления, компьютерное моделирование. 


Введение. Для исследования летательных аппаратов (ЛА) сложной геометрической конструкции 
необходимо создавать математические модели (ММ). Обычно они основываются на разработан- 
ной Эйлером ММ свободного движения твёрдого тела. Основная трудность её идентификации со- 
стоит в определении действующих на ЛА сил (лобовой, подъёмной и боковой), а также возника- 
ющих при этом моментов вращения. Для каждого конкретного ЛА эта задача может быть решена 
только экспериментально — с использованием аэродинамической трубы или компьютерного мо- 
делирования. 

Авторы принимали участие в разработке транспортно-коммуникационной системы МААТ (в 
Седьмой рамочной программе, финансируемой Евросоюзом). Следует отметить, что проект МААТ 
рассматривает систему двух типов аэростатных ЛА: крейсера, постоянно находящегося на страто- 
сферных высотах, и челноков, доставляющих грузы и пассажиров на крейсер и обратно. Послед- 
ние работают в наиболее сложных условиях, связанных с чрезвычайно широким диапазоном из- 
менения атмосферных параметров. Таким образом, авторы столкнулись с необходимостью по- 
строения математической модели транспортного челнока — дирижабля, курсирующего между 
аэропортом и дирижаблем-крейсером, барражирующим в стратосфере на высоте 15—20 км. В дан- 
ной статье описывается весьма сложный и ответственный этап исследования — идентификация 
силовых и моментных характеристик одной из конструкций челнока. Причём, как принято при си- 
стемном решении научных задач, на первом этапе целесообразно запланировать изучение иде- 
ального варианта подъёма челнока — всплывание за счёт архимедовых сил с доступными при та- 
ком принципе движения скоростями. Необходимо учесть также, что в реальных условиях ориен- 
тация челнока может меняться под воздействием внешних возмущений. 

Описание конструкции исследуемого варианта челнока. Модель челнока изометрически 
изображена на Рис. 1. Для разработки её 3З)-чертежа и математического описания используется 


1 Работа выполнена в рамках инициативной НИР. 
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принятая в теории воздухоплавания ортогональная система координат. Ось ОХ — в горизонталь- 
ной плоскости и по направлению основного горизонтального движения. Вертикальная ось ОУ — 
вертикальное направление. Рассмотрение боковых движений челнока привязывается ко второй 
горизонтальной оси — ОХ. 





>“ мрй < 


Рис. 1. Модель челнока системы МААТ 


Челнок в исследуемом варианте представлен жёсткой каркасной конструкцией круглого 
в плоскости ОХС сечения. В верхней части конструкции под куполом по замыслу проектировщиков 
помещена камера переменного объёма с эластичными упругими стенками. На максимальной высо- 
те плавания она заполняет весь объём корпуса, кроме помещённой внизу гондолы-кабины (также 
круглого горизонтального сечения). Более подробно конструктивные элементы челнока не прора- 
батывались, т. к. предполагались его аэродинамические исследования — в первую очередь с целью 
определения размещения и выбора мощностных характеристик двигателей. 

Физические и конструктивные характеристики модели приводятся в табл. 1. Некоторые 
из этих параметров носят прикидочный характер, т. к. разработка представляла собой эскизный 
проект. 

Таблица 1 


Физические и конструктивные параметры челнока 


_ Диаметр, м Положение центра масс и плавания 
Параметр | Масса, кг | Объём, м? | Высота, м р Р 


25000 583783 88,78 128,11 78,62 -17,97/-63,54 


Постановка задачи. Таким образом, на основании данных о предполагаемой конструкции чел- 
нока и о среде его плавания во всём диапазоне необходимо оценить значения возникающих при 
движении сил и моментов, требующих активного вмешательства со стороны систем управления. В 
качестве активно задаваемых входных параметров эксперимента выбраны три наиболее важных 
по предварительно сделанным оценкам: 

1) высота местонахождения челнока над уровнем моря; 

2) скорость вертикального подъёма челнока; 

3) угол наклона челнока относительно горизонта (в силу симметричности форм и верти- 
кальности подъёма теряют смысл привычные для авиаторов понятия крена и дифферента — они 
равнозначны в этом случае). 
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За диапазон пространственного исследования принят высотный диапазон от точки старта 
(предположительно на уровне моря) до точки стыковки со стратосферным крейсером (предполо- 
жительно на уровне 15 км). Предварительные исследования требуемого и возможного быстродей- 
ствия челнока показали, что диапазон возможных скоростей его всплывания ограничен значени- 
ем 7-8 м/с. Наконец, диапазон возможных углов наклона челнока для полноты картины исследу- 
емого феномена целесообразно взять несколько меньше максимального — 90°. Последний, разу- 
меется, в реальном полёте немыслим, т. к. находится на границе начала полного переворота ап- 
парата — оверкиля. Однако свойства исследуемой конструкции в экстремальных режимах, при- 
ближённых к критическому, могут быть использованы при её совершенствовании. 

В соответствии с заявленной целью исследования необходимо, варьируя выбранные вход- 
ные переменные (их диапазон оговорён при постановке задачи), получить математическое описа- 
ние сил и моментов, возникающих при всплывании челнока. При исследовании выбранного иде- 
ального движения отмечено три вида силовых воздействий на аппарат: 

1) сила сопротивления всплыванию (антипод подъёмной силы); 

2) сила бокового смещения, возникающая при наклонном положении челнока из-за его 
взаимодействия с набегающим при всплывании потоком воздуха; 

3) вращающий (опрокидывающий или возвращающий к исходному состоянию) момент, 
формируемый взаимодействием аппарата с набегающим при всплывании потоком воздуха. 

Иными словами, необходимо построить математическую модель аэродинамических 
свойств исследуемого варианта челнока. 

Предварительное планирование эксперимента по исследованию характера силовых и 
моментных воздействий на челнок среды плавания. В соответствии с поставленной зада- 
чей разработан план эксперимента, числовые характеристики которого представлены в табл. 2. 
Таблица 2 
Числовые характеристики плана эксперимента 


Параметры всплывающего челнока 
Параметры планируемого эксперимента 
Угол наклона челнока, град. | Скорость всплывания, м/с | Высота плавания, км 


Нижний уровень ЗН ЗА ОЧ ЗАЕЕЛО ОНИС ОИ 
Диапазон варьирования 15-1.0-10=55 





Верхний уровень ОС И ПО И ООО 


Сравнительно грубая шкала исследования объясняется высокой трудоёмкостью процесса 
электронной продувки с помощью компьютерного моделирования и прикидочным характером ис- 
следования. Переменный шаг квантования уровня изменения скорости связан с тем, что предва- 
рительно в проекте выбрано её ориентировочное значение — 5 м/с. 

Полученная в результате сетка 6х5х4 определяет количество опытов — 120. 

СотриваНопа/ ЙА дупат/с$ (СЕО) — это специальное программное обеспечение (СПО), 

моделирующее течения жидкостей и газов с учётом турбулентности, теплообмена и химических 
реакций. Применение данных программ позволило получить искомые значения сил и моментов. 
При этом использовалось свободно распространяемое ПО, доступное на платформе [лих. Были 
проведены опыты по компьютерному моделированию в среде СПО, позволившие получить чис- 
ленные оценки результатов запланированных опытов. 
Математический аппарат исследования силовых и моментных воздействий на челнок 
среды плавания. Аппроксимация полученных экспериментальных зависимостей осуществляется 
обычно с использованием аппарата регрессионного анализа. Структура искомого уравнения ре- 
грессии (УР) принимается в подавляющем большинстве случаев как традиционная полиномиаль- 
ная форма, общий вид которой задаётся следующим выражением: 
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П П П 
у(х) => В -х-+ УВ, -х,-х,+ >» Ви -Х,-Х, -Хич... (1) 
7—0 7-0 ИрК=о 
Здесь у(х) — выходная величина искомой зависимости, т. е. в нашем случае силы или моменты 
Х! 
воздействия среды плавания на челнок; х =| х, | — вектор входных переменных, т. е., согласно 
х 


3 
табл. 2, угол наклона челнока — х, , высота его плавания — х, и скорость подъёма — х,. Пара- 


метрами УР являются константы БВ, 


я назначение которых состоит в обеспечении максимальной 


адекватности выражения (1) экспериментальным данным. При этом Б, — коэффициенты линей- 
ной части УР; х, — единичный вектор размерности № =120 (количество опытов); В, — свобод- 


ный член УР. Во второй сумме (1) 6, — коэффициенты членов УР второго порядка при парных 


произведениях переменных х,-х, (при / = 7 формируются квадраты переменных). Третья сумма 


формирует нелинейные члены УР третьего порядка по такому же принципу, ит. д. 

Конкретный порядок описания выбирается индивидуально для каждого исследуемого яв- 
ления или объекта, т. к. регрессионное уравнение типа (1) не отражает фактической физической 
структуры взаимодействия параметров, определяющих процесс формирования аэродинамических 
сил. Поэтому его компоненты и их коэффициенты формируют лишь количественную составляю- 
щую исследуемого явления, что часто приводит к невысокой точности математического описания. 
Коррекция структуры и параметров эксперимента по исследованию свойств челнока. В 
связи с вышеизложенным важно проанализировать качество полученного математического опи- 
сания и на его основе критически оценивать запланированный эксперимент, выявлять парамет- 
рические причины невысокой точности модели (1), меняя диапазоны варьирования выбранных 
параметров либо даже сами параметры. 

Так, при решении поставленной в работе задачи на этапе регрессионного описания иско- 
мой зависимости выяснилось, что модель (1) даёт очень низкую точность описания при опти- 
мальных по МНК (методу наименьших квадратов) коэффициентах регрессии. Отклонения в ряде 
точек превышали 100 %, что недопустимо даже при прикидочных инженерных расчётах. 

Математический анализ причин данного феномена проведён путём включения — исклю- 
чения нелинейных компонентов выражения (1) и оценки их влияния на точность математического 
описания силы сопротивления всплыванию челнока. Анализ выявил два фактора формирования 
погрешностей: 

1) излишне широкий диапазон варьирования угла наклона (большинство ошибок и наибо- 
лее значительные из них приходились на опыты с углом 90°); 

2) плохое качество описания зависимости силы всплывания от высоты. 

Последнее связано с экспоненциальной зависимостью от высоты плотности воздуха, силь- 
но влияющей на аэродинамику. Полином плохо аппроксимирует экспоненту. 


Таблица 3 
Скорректированные числовые характеристики плана эксперимента 


Параметры преобразованного Скорректированные параметры всплывающего челнока 
эксперимента 


75 
жний уровень О 15 
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Поэтому условия эксперимента были скорректированы (табл. 3). Исключён уровень варьи- 
рования угла 90°. Параметр «высота» заменён параметром «плотность» в функции высоты, соот- 
ветствующей принятым её уровням. 

Таблица 4 
Результаты имитационного моделирования 


Параметры опытов 
челнока, а ха, Ре вания, а Ре .Н чение, Н 


а 
О т 
м | 1500 | 125 | 250 | 390 | 39%8я | 007 — 
2 | 1500 | 1250 | 40% | ивыво | снизу | 590462 
| 1500 | 1250 | 50% | 530 | 970562 | 0009 — 
| 1500 | 1250 | 60% | 24840 | 543565 | 30002 


В 
| 3000 | 045 | 7500 | ывмо | 558069 | 5905976 


И о Ш 
в | 6000 | 045 | 500 | мо | 56862 | 0460 — 
| 6000 | 05 | 60% | 29м%р | 526605 | 01089 — 
ЕЯ ПАВ ЖИ ПОЕТ: ПИТЬ 


Результаты эксперимента и их исследование. В соответствии с планом (см. табл. 3) прове- 
дено имитационное компьютерное моделирование процесса всплывания челнока в идеально не- 
подвижной среде. Результаты сведены в табл. 4. Здесь представлены фрагменты полных данных 
(всего проведено 120 опытов). В третьей колонке — значение плотности воздуха на высоте, соот- 
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ветствующей уровню варьирования (табл. 2). В пятой колонке приводится результат измерения 
силы всплывания. 

Поскольку выходная переменная рассматривается как функция трёх независимых пере- 
менных, регрессионный анализ проводился относительно полинома третьей степени. При таком 
порядке, как известно из [2], регрессионная полиномиальная аппроксимация начинает проявлять 
наряду с фильтрующими свойствами свойства генератора шума. Этот эффект проявляется следу- 
ющим образом. На фоне увеличения точности аппроксимации опытных точек и общего уменьше- 
ния среднеквадратической ошибки (СКО) возникают значительные отклонения аппроксимирую- 
щей функции от огибающей этих точек. Это приводит к ухудшению качества аппроксимации. По- 
этому полученный исходный полный полином третьего порядка, представленный в (1) первыми 
тремя суммами, был исследован с точки зрения характера влияния всех составляющих на струк- 
туру и точность. Незначимые и искажающие результат члены отбрасывались. 

Результат обработки уравнения (1) представлен в табл. 5. Из 20 членов полного полинома 
3-й степени при трёх независимых переменных в УР оставлено 14. Их структура (определяется по 
индексам обозначения) и численные значения коэффициентов сведены в табл. 5. 

Таблица 5 
Коэффициенты УР 
р р 





Результаты расчёта силы сопротивления всплыванию по сформированному УР и относи- 
тельной ошибки оценки этой силы приведены в табл. 4 (шестая и седьмая колонки соответственно). 
При такой структуре УР дало наименьшую ошибку аппроксимации по сравнению с осталь- 
ными вариантами: дисперсия относительной ошибки составила 0,0198, а среднеквадратичное от- 
клонение — 0,1408. Анализ характера ошибок показал, что они имеют нормальное распределение 
вероятностей. Поэтому доверительный интервал для 67 % вероятности невыхода ошибки состав- 
ляет ^0,14, а для 95 % -0,28. 
Само УР в принятых выше обозначениях имеет вид 
у (х) =-19322,7 + 419,75. х, + 23498,3.х, +14965,7.х, -399,29.х, .х, — 
—26508,5.х,.х, -1653,27. х; +547,71.х,.х,.х. -505,9.х, -Х: + (2) 
—35,85. хх: -4,05.х2.х, - 3416,7. х, Хх; +1,07.х2.х, +12601,6.х$. 
Авторы исследовали специализированный летательный аппарат аэростатного принципа всплыва- 
ния. На рис. 2 построены графические З)-иллюстрации зависимости силы сопротивления его 
всплыванию от скорости и угла поворота для различных высот. 
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Рис. 2. Зависимость силы сопротивления воздуха вертикальному всплыванию челнока МААТ на высоте: а — 0 км, 
6 — 5 км, в— 10 кмиг— 15 км 
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Выводы. Полученные результаты демонстрируют эффективность имитационного компьютерного 
моделирования и регрессионного анализа при построении математических моделей летательных 
аппаратов. В ходе реализации проекта МААТ эти результаты будут использованы для разработки 
систем автоматического управления полётом. Невысокая точность математической модели силы 
всплывания компенсируется в последующих исследованиях совершенствованием специализиро- 
ванных программ. 
Библиографический список 

1. Гтрасе ог {Пе Геедег Аегодупаптс$ СПагасеп$с$ оп {Пе Ромег оГ Сотго! Асйоп$ п Зеаау 
апа Тгапчепе Кедитез / \. РУМКПором, [её а|.] // Аегозрасе Еесгогс апа Амотс$ бу$етз Сопегепсе 
(Рпоетмх, Айтопа, ЦА, Ос. 30 — Мом. 1, 2012). — Ма. 3. — Р. 140-150. 

2. Юл, Дж.-Э. Теория статистики / Дж.-Э. Юл, М.-Дж. Кендэл ; пер. с англ. Ф. Д. Лившиц [и 
др.]. — 14-е изд., пересмотр. и расш. — Москва : Госстатиздат ЦСУ СССР, 1960. — 779 с. 

3. Смирнов, Н. В. Курс теории вероятностей и математической статистики для технических 
приложений / Н. В. Смирнов, И. В. Дунин-Барковский. — 2-е изд., испр. и доп. — Москва : Наука, 
1965. — 511 с. 


Материал поступил в редакцию 07.03.13. 


ВеГегепсе$ 

1. М. Р5МКПором, ЕЁ а|. Ттрасе ог {Пе Геедег Аегодупапт!с$ СПагацет$йс$ оп {Пе Ромег о Соп- 
{го| АсНоп$ п Зеау ап Тгапзепе Кедитез. Аегозрасе Еесгопс апа Амогс$ бу%етз Сопегепсе 
(Рпоетмх, Айтопа, Ц$А, ОсЕ. 30 — Мом. 1, 2012), ма. 3, рр. 140-150. 

2. УшЕе, С. Ч., Кепдай, М. С. ТПе ТВеогу ог За!$8сс ; {гап$1. Гот ЕпдИ$й Бу ЕО. ИМм$ПК$ её а1. 
1 А гембеа апа епагдеа ед. Мо$соми : бСоз{ай2аае С25Ц $$$, 1960, 779 р. (т Кизчап). 

3. типом, М. М, Бипт-ВагКоу$Ку, Т. М. Киг$ {ео уегоуапо$е] {1 таетайспезКо] Уай$ик 
Чуа фехиспезЮКх ри о7Пепп. [ТПеогу ог сПпапсез ап та{Петайса! $аИ$Нс$ соигзе Гог епдтееппда ар- 
рИсаНоп$.] 2п4 гембе4 апа етагде4 ед. Мо$соми : МацКа, 1965, 511 р. (т Кизчап). 


КЕЗЕАВСН ОМ ВЦОУАМТ ЕОВСЕ ОЕРЕМОЕМСЕ ОЕ $РЕСТАТТ7ЕО АЕВОЗТАТ ОМ 1Т$ 
МОТТОМ УТВАВЕЕ$: 


В. А. М№еудом, У. 1. $олАа 
(Роп Чае Тесптса! Упмег®Ку) 


Тре Ёегпанопа!/ зсепИИс-гезеагсй МААТ (Ми/ЕБоду Адуапсеа Ар Гог Тгапзрой.) ргодесЕ зирропеа Бу {пе Еи- 
гореап Ипюп оп {пе $еуеп т Гатеиогк ргодгат зо№ез5 {пе д/оБа! рго Мет оп {Те есоодгса/У 5за!е гапзропбайоп. 
Опе о! [5 есноп$ Ваз [Ве айт Ю гезеагср те аегодупат/!с$ о! Те 5реса/геа ай-аиае БаШооп$ — хо-са!еа 
5ПШИШЕ5. ТПРеу аге аез!дтеа Гог 5иррогипа соттигсайоп$ Бемееп {Ве еапйП-Базеа а!ро!5 апа {Ве аего&аЁ5 
и//СВ аге регтапепйу сги/та т {Ве гайюзрпеге — сгибег5. Опе о! те аегодупат/с гежеагсВ {а$К$ пеедеа ю 
деуе/ор ап арргорпае та етаЕса!/ тосе! 15 5о/уед. Тре ргоМет оп Те таетайса! ЧезсирИоп о! {пе БиоуапЕ 
Гогсе спапде$ ипаег {Пе пи Ше р/ипата т #5 Бо ипазигЬеа апа @бигреа ро$!оп [5 $ЕЁ апа зо/ед. Гог {пав 
а хтиавоп ехрейтепЕ Базеа оп {Пе хо-са!еа “еесгоп/с айПои/” {пгоидВ те сотршег тодетод 1/5 сопаисеа. 
Тре сотрийпд гезий оЁ те Биоутд агад Гогсе ипдЕеГ уагутда сопа!Ноп$ аге ргосез$е$ пгоидй {пе гедге$юоп 
апа/у5/5 &есптаиев. Тре оатеа ехргез$/оп рептт!&5 {о аепёу {пе та етайса!/ тодЕ! о! {пе эпие тобоп ду- 
пат/с$ ргезещеа {пгоидв Ешег$ едиайоп$. ТПе гези!5 аге из гайеа Бу те вЫ апа дгар!Р/са! аеа. 

Кеуигога$: раюзрПЕге, а/г ие СЕ, аего%аё Биоута рипс/рИЕ, геячапсе Гогсез апа тотепё, сотриг ути/авоп. 


1 Тре гезеагсН {5 допе минт {Пе гате оЁ {Не тдерепаепе В&р. 
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